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ON THE CALCULATION OF FLOW ABOUT OBJECTS 
TRAVELING AT HIGH SUPERSONIC SPEEDS 
By A. J. Eggers, Jr. 


SUMMARY 


A. procedure for calculating three-dimensional steady and nonsteady 
supersonic flows with the method of characteristics is developed and dis- 
cussed. The flow is assumed to be adiabatic and inviscid, although it 
may be rotational and the gas may exhibit both thermal and caloric imper- 
fections. The latter features of generality are retained in the analysis 
since it is known that the phenomena associated therewith may signifi- 
cantly inf luence flows at high supersonic airspeeds . A further study of 
the compatibility equations holding along characteristic lines reveals 
that at Mach numbers sufficiently large compared to 1, flow in the oscu- 
lating planes of stre aml ines may, in regions free of shock waves, often 
be of the generalized Prandtl-Meyer type. Surface streamlines may, under 
such circumstances, be approximated by geodesics. These results hold for 
nonsteady as well as steady flow, provided only slender shapes are con- 
sidered, and provided the induced curvature of the flow associated with 
the nonsteady motions does not exceed in order of magnitude the total 
curvature of the flow. Steady two-dimensional-flow equations may thus 
be applicable to a wider class of flows, and hence shapes, at high super- 
sonic speeds than was heretofore thought. 


INTRODUCTION 


The calculation of flows about objects, primarily missiles, travel- 
ing at high supersonic speeds is now generally accepted as a matter of 
more than academic interest. The difficulty of these calculations stems 
in large part from the fact that at such high speeds disturbance veloc- 
ities are not necessarily small compared to the velocity of sound, nor 
are entropy gradients necessarily negligible in the disturbed flow field 
about a body, even though it may be of normal slenderness. Thus, for 
example, the relatively simple linear theory, which haB proven so valu- 
able in studying flows at low supersonic speeds, loses much of its utility 
in the study of high-supersonic-speed flows. In the quest for methods 
especially stilted to calculating high-supersonic-speed flows, notable 
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progress has been made in the development of similarity laws relating 
the flows about slender three-dimensional shapes in both steady (see 
references 1, 2 , and 3) and nonsteady motion (see references 4 and 5)* 
Steady two-dimensional flows have received perhaps the greatest attention 
from the standpoint of calculating specific flow fields, and it would 
seem that with tools ranging from the characteristics method (see, e.g., 
references 6 and 7) to the generalized shock-expansion method (refer- 
ence 7) the problem is reasonably well in hand, at least insofar as 
inviscid, continuum flow is concerned. A more or less analogous situa- 
tion exists with regard to the nonlifting body of revolution (see, e.g., 
references 6, 8, 9, and 10) although it seems that only in the case of 
the cone has a method (reference 10) of simplicity comparable to that of 
the linear theory been developed for calculating the whole flow field. 

When one departs from these relatively simple flows, the number of • 
tools for carrying out practical calculations decreases sharply. Thus, 
for example, in the category of inclined bodies of revolution, it appears 
that only bodies at small angles of attack have been handled adequately, 
usually by either the method of characteristics or seme other step-by-step 
calculative procedure (see, e.g., references 6, and 11 through 14). In 
the case of steady flow about general three-dimensional shapes, aside 
from Newtonian flow concepts, which are strictly applicable at Mach num- 
bers exceeding all limits, only the characteristics method has apparently 
thus far received serious attention (references 15 > 16, 17, and 18). 

Sauer's treatment of this method (reference 18) is especially neat, 
entailing only the assumption of ideal (i.e., thermally and calorically 
perfect) gas flow and yielding compatibility equations in a relatively 
simple form. Application of the method, although it would undoubtedly 
prove tedious and time consuming, is, as pointed out by Sauer, certainly 
feasible with the aid of present day high-speed computing machines. It 
is clear, of course, that the relatively exact solutions obtainable with 
the method of characteristics provide an invaluable check against the pre- 
dictions of more approximate but perhaps simpler theories. Indeed, as 
demonstrated in reference 10, a study of the compatibility equations of 
the characteristics method can prove useful in determining simplified 
methods for calcinating more complex flow fields. 

With these points in mind, it is undertaken as the first principal 
objective of the present report to reconsider the characteristics method, 
particular attention being given to its application to high- supersonic - 
speed flows in which, aB is often the case, air does not behave as an 
ideal gas. The assumption, then, of ideal gas flow is relaxed, and to 
preserve insofar as is possible the element of simplicity in the method, 
pressure and flow inclination angles are employed as dependent variables 
(see reference 7) rather than, for example, velocity components as were 
used by Sauer. Extension of the method to the study of nonsteady flows 
is also considered; however, the remaining principal objective of this 
paper is to show how results of the characteristics theory can be exploited 
to deduce a simplified procedure for calculating high-supersonic-speed 
flows of both the steady and nonsteady types. 
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ANALYSIS AND DISCUSSION 


Compatibility relations describing the behavior of fluid properties 
along characteristic lines in supersonic flow may, of course, be obtained 
by proceeding formally with the theory of characteristics for the quasi- 
3 in ear partial differential equations which depict the flow. In the 
interests of simplifying both the derivation of these relations and their 
resultant forms, however, it seems desirable to proceed in a more intuitive 
manner, assuming a priori that the characteristic lines are Mach lines and 
streamlines (in steady flow), and utilizing the implication from two- 
dimensional- flow studies that perhaps the most convenient dependent vari- 
ables are pressure and flow inclination angles. With this approach in 
mind, we employ the Euler momentum equations. 


gU lr |2 + v |SU lf |!U ,i|£ 

at ox Sy Sz P ox 

(1) 

51 + u M + v Sv + w Sv , . 1 

St Sx Sy Sz P Sy 

(2) 

St Sx Sy Sz P Sz 

(3) 

the continuity equation. 

Sp , S(pU) . S(pV) S(pW) _ A 
St Sx Sy Sz 


the equation of state, 

P = p(Pjs) 

\Sl 

and the energy equation. 

|2. + v |£. + y |S. + w is. . 0 

St Sx Sy Sz 

(6) 


where U, V, and W are the components of velocity at time t along 
the x, y, and z axes, respectively, of an element of the fluid of 
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density p, static pressure p, and entropy s 1 (see appendix for list 
of symbols). To put these expressions in a more tractable form, it is 
convenient to aline the x axis at time t with the direction of the 
resultant velocity at the origin of the coordinate system. Thus equa- 
tions (l) through ( 4 ) and equation ( 6 ) simplify, respectively, in the 
region of the origin, to 


^ . TT 5U . 1 . 

r — + U — + — tt*- - 

at ax P ax 
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and 
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+ = 0 
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( 11 ) 


’ which relations are basic to the 


subsequent analysis. 


Steady Flow 


Characteristics method .- It is clear that in this case all derivatives 
with respect to time disappear from the above relations. Thus, assuming 
there are no shock waves present in the region of the origin, 2 we may 
write, with the aid of equations (5) and (11), 


1 For certain calculations it may be desirable to proceed from more gen- 
eral equations which include effects of heat and mass addition to (or 
subtration from) the flow as well as effects of impressed forces (e.g., 
gravitational or magnetic). Such a procedure may easily be developed 
from that presented here by following the method of Guderley (refer- 
ence 6) for two-dimensional flow. 

2 If shock waves are present, the appropriate oblique shock equations are 
employed. 
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dp _ dp 

dx dp 

■where a is the local speed of sound in the fluid. Combining 
equations (7 ), (10) , and (12), there is then obtained the relation 


k = -L£R (12) 

dx a 2 dx 


dp _ -pU g f* 1 /dV + dW V 
dx M 2 -l _ U \ dy dz ) _ 


(13) 


or, defining A as the angle between the x axis and the tangent to the 
projection of a stre aml ine in the x - y plane, and, in an analogous 
manner, the angle 5 in the x - z plane (see sketch), we have 


Streamline 


2 



dp _ -pU 2 f oA [ d.S N 
dx M 2 -l \ dy dz ) 




Transforming the derivatives with respect to x and z to derivatives in 
the characteristic or C lz and C 2Z directions in the x - z plane (Ci z 

is positively inclined with respect to x, thus d()/dx = [M/(2 a/M 2 -1)] 
[d()/dC lz + d( ) /dCgg] and d()/dz = (M/2)[d( )/dC lz - d()/dc 2z ] ), there 

results from this equation 


dp 

+ dp _ 

PU 2 

d6 

dS 
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dCiz 

dCsz 
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’ dc lz 
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(15) 


In an analogous manner, there is obtained from equation (9) the relation 

(16) 


dCxz 


- dp = -PU 2 ( 

dCs z */M 2 -l \ 


d6 dS 
dCj_ z dc^j 
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Adding these two expressions then yields 


8Ci Z 

while subtracting yields 


85 

^ c iz 


— f'— ^ 

M V3 y / _ 


u I 2 r 8s _ 1 Y 8 a \ ‘ 
[Cl L 8c PZ M \ 8y / _ 


Equations (17) and (l8) are compatibility equations for charateristic 
or Mach lines in the x - z plane. 3 Indeed, if it is further required 
that the x - z plane be the osculating plane of the streamline passing 
through the origin, that is, the plane containing the principal radius 
of curvature and tangent to this streamline (at the origin), then these 
equations are the essential relations for deter mining pressure and flow 
Inclination throughout a flow field. This point becomes evident when 
it is observed that with the imposed requirement (viz., 8a/8x = 0), the 
additional information derived from studying flow in the x - y plane 
is simply that deduced from equation (8), or, as would be expected, 

|| = 0 (19) 

In proceeding to construct a flow field, however, it is necessary to 
know the manner in which the osculating plane rotates and, correspondingly, 
how the principal curvature varies along a streamline. This information 
is easily obtained from equations (2) and (3) by differentiating with 
respect to x, thus yielding 


I 8 /SjjjA 3 a 85 

pU 2 8x \ 8y/ 3z 3x 


1 _ 3 _ 
pU 2 3x 



\ 3z u 


respectively. 

With the aid of these and the previously derived expressions, we now 
consider how the characteristics method for steady three-dimensional flows 


3 It is noted that these expressions contain not only derivatives in the 
characteristic directions but also derivatives with re*spect to the 
independent variable y. This type of result is to be expected as 
ointed out by Coburn. (reference 19). 
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can be applied. In particular, only the initial value problem will be 
considered here. Thus it is first assumed that flow conditions are 
known, or in some manne r may be determined on a surface 4 in the disturbed 
flow field to be investigated. We then consider flow in the region of 
the osculating plane of one of the streamlines passing through this sur- 
face (see sketch below). 



Using known conditions at points A and B, pressure and flow inclination 
at C are calculated employing equations (17) and (l8) as difference 
equations. Orientation of the osculating plane and curvature of the 
streamline at C are then determined with the aid of equations (20 ) 
and (21) and the knowledge of x and 86/dx. Other properties, Buch as 
velocity, temperature, and Mach number are determined at point C in a 
manner analogous to that for two-dimensional flow (see, e.g., refer- 
ence 7)* Once these calculations have been carried out at a number of 
points, like C, close to the initial surface, then the whole procedure 
is repeated, progressing from these new points (or points interpolated 
therefrom) of known flow conditions to other points farther removed from 
the surface. Thus the flow field is constructed moving downstream (or 
upstream) from the initial surface. 

Variations on these calculations are frequently required in practical 
applications; for example. If a shock wave is encountered in the construc- 
tion; it is usually necessary to solve simultaneously the equations for 
the shock and the compatibility equation for the characteristic line 
intersecting the shock in order to determine flow conditions at a point 
adjacent to the shock. With an analogous procedure, points on the sur- 
face of a body may be treated. In any case, it would appear that 


4 If no other Information on the flow field Is available, this surface 
cannot, as is known from the theory of characteristics, be a character- 
istic surface. However, as is frequently the case in aerodynamics, this 
restriction is eliminated since the initial surface intersects a shock 
wave or body of known shape, or both, which result provides additional 
boundary conditions at the terminal edges of the surface. 
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complication of calculation has been reduced by working with pressure 
and flow inclination angles as primary dependent variables, no loss in 
generality being sustained insofar as restrictions on the thermodynamic 
behavior of the gas is concerned. The gas may exhibit both thermal and 
caloric Imperfections without invalidating the previously developed 
equations, the extent to which these imperfections are manifest influenc- 
ing only the form of the equation of state of the fluid and the relations 
defining its specific heats. 5 The provision for caloric imperfections 
would appear (see reference 7) to be especially desirable when using the 
method to calculate hypersonic flows. 

Application of the method to the calculation of a specific flow 
field hinges on the determination of flow conditions along the initial 
surface . This is a separate problem, the solution of which has thus far 
at least been special to the particular type- of flow under consideration.® 
Indeed, only the airfoil has apparently been treated rigorously In this 
connection without the restriction that air behaves as an ideal gas (see 
reference 22). 

Imperfect gas flow downstream of the throat of a hypersonic nozzle 
could, of course, be analyzed by the characteristics method of this 
paper. With the restriction to Ideal gas flows, an additional application 
of importance also suggests itself, namely, to the calculation of flows 
about inclined bodies of revolution at angle of attack (not necessarily 
small). In this case flow conditions along a surface close to the vertex 
of the body can be calculated approximately with the aid of references 25 
and 26, the accuracy of the approximation increasing with the closeness 
of the surface to the vertex. Flow downstream of this surface can then 
be determined after the manner discussed, using equations (17), (l8), 

(20), and (21) in the reduced forms 


Jil. 
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6a 5s 
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TpM 2 8x \ 6; 

T' 
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5 Since only flows of dense air are treated here, heat capacity lag 
effects are considered negligible. 

e In this connection see, for example, the work of Crocco (reference 20), 
Munk and Prim (reference 21), and Kraus (reference 22) on the two- 
dimensional airfoil problem, and Shen and tin (reference 23) and 

Cabannes (reference 24) on axially symmetric flow about bodies of 
revolution. 
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» and 

v d 2 S _ l 8 / 8p \ [~ 8s 1 8p (. jz __ » 1 85 

8 x 2 7 pM 2 3x \ 8 z / Uz 7 pM 2 8 x 8 x 

respectively, where 7 , the ratio of specific heat at constant pressure 
to specific heat at constant volume, would be considered const ant (equal 
to approximately 1.4 for air). 

Let us now turn our attention to the consideration of a more 
approximate, and by the same token more simplified, method of calculat- 
ing the steady three-dimensional flow of a gas at high supersonic speeds. 

Approximate method .- It is well at the outset of this analysis to 
establish, insofar as is practicable, the type of flows to be treated. 

In this connection, it is convenient to employ the hypersonic s imil arity 
parameter (i.e., the product of the flight Mach number and say the thick- 
ness ratio of a body) as a measuring stick. In flows characterized by 
values of the hypersonic similarity parameter small compared to 1 , that 
» is, flows in which the body is extremely slender and lies close to the 

axis of the Mach cone, there is no apparent reason to believe that the 
linear theory will not be as useful an approximate method of calculation 
' as at low supersonic speeds. In flows characterized by values of the 

parameter up to about 1 , the second-order theory first enunciated by 
Busemann (reference 27) for airfoils and more recently generalized to 
three-dimensional flows by Van Dyke (reference 9) and Moore (reference 28) 
should prove a useful approximation. On the other hand, for flows about 
more or less arbitrary shapes, there is apparently no approximate method 
of calculation generally applicable with engineering accuracy at values 
of the hypersonic similarity parameter appreciably greater than 1 . 

In the limiting case of indefinitely high free-stream Mach number 
(and hence similarity parameter) and a ratio of specific heats equal to 1 , 
we have the Newtonian impact theory (reference 29) and its refined 
counterpart, accounting for centrifugal forces in the disturbed flow, 
developed first by Busemann (reference 30) and more recently treated 
by Ivey, Kluhker, and Bowen (reference 31) • The impact theory has been 
employed with some success by Grimminger, Williams, and Young (refer- 
ence 32 ) and others to predict surface pressures on bodies of revolution 
at values of the similarity parameter appreciably greater than 1 , although 
it should be remarked in passing that this success is in part, at least, 
fortuitous, as perhaps is best evidenced by the fact that the more exact 
theory (within the framework of the underlying assumption of M— >«, 7 — >1) 
of Busemann is considerably less accurate under corresponding circumstances. 
" As Bhown in reference J } neither the Newtonian impact nor the Bus emann - 

theory apply with good accuracy to airfoils except at values of the B ind - 
larity parameter quite large compared to 1 , correspon din g, for example, 
in the case of thin airfoils to flight speeds considerably in excess of 
the escape speed at sea level. Perhaps the foremost shortcomi ng of these 
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theories is, however, that, irrespective of . hhe shape to which they are 
applied, they provide no information on the structure 7 of the disturbed 
flow field which is, of course, of finite extent adjacent to the surface 
at flight Mach n umb ers presently of interest (say Mach numbers less than 
the escape Mach number at sea level). Such information is, for example, 
important to the determination of the flow about control surfaces and 
the like which may be located in this field. 

In view of the preceding discussion, it seems clear that in the high- 
supersonic-speed flight regime, the need for an approximate method of 
anal ysis lies in the realm of flows characterized by values of the hyper- 
sonic s imi larity parameter greater than 1. An attempt will therefore be 
made to obtain a method meeting part of this need, attention being focused 
primarily on flows characterized by large values of the similarity param- 
eter. To this end, it is convenient first to employ equation (l4) rewrit- 
ten in the form 


where 




D - d&/dCiz 

z Ss/dc 2z 


( 22 ) 


(23) 


Now consider for the moment a surface streamline alined in the x direc- 
tion, and impose the requirement that the x - z plane be tangent to 
this streamline and normal to the surface at the point of tangency (the 
origin). The x - y plane is then, of course, tangent to the surface 
at this point. Observing the last term in the brackets on the right-hand 
side of equation ( 22 ), it is noted (see sketch) that 



7 This consequence is traceable primarily to the assumption of 7=1 
which leads to the well-known result that the disturbed flow field is 
confined to an infinitesimal region adjacent to the surface of a body. 
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where r is the radius of curvature of the line normal to the projec- 
tions of streamlines in the x - y plane, and passing through the origin. 
At the high Mach numbers under consideration, the disturbed flow field is 
confined to a region of small extent normal to the surface of a body; 
hence it may be expected that r will be primarily a function of body 
shape and atti tude. 8 This being the case, it follows then that the 
term (l/ */m 2 -1) (l/r) will decrease in absolute magnitude with increasing 
Mach number of the flow about the body. Consider now the term 
(d8/dx)(l-D z )/(l+D z ). We note that 3&/8x = l/R where R is the radius 
of curvature of the projection of a streamline in the x - z plane and, 
by reasoning analogous to that used in considering r, is not expected 
to vary significantly with Mach number in the disturbed flow field. Let 
us assume for the moment that the quantity (l-D z )/(l+D z ) is also rela- 
tively independent of Mach number. With this assumption, it is clear 
that equation (22) approaches the equation for two-dimensional flew as 
the free-stream Mach number, and hence the hypersonic similarity param- 
eter of the flow becomes large compared to 1. The compatibility equations 
(equations (17) and (l8)) are affected in a similar manner; thus it is 
apparent that the flow when viewed in the x - z plane approaches the 
two-dimensional type. In this case, however, as Bhown in reference 7* 
so long as the Mach number and ratio of specific heats of the disturbed 
fluid are not too close to 1, D z is small compared to 1, and hence the 
flow approaches the generalized Prandtl-Meyer type (i.e., flow in which 
pressure and inclination angle are approximately constant along curved 
first-family Mach lines ) . Our flow equation may then be written 



where it is required explicitly that 

85 >;> 1 

dx a/M 2 -1 dy 

or, in effect, that disturbances associated with the divergence of stream- 
lines in tangent planes must be of secondary importance compared to those 
associated with the curvature of streamlines in planes normal to the sur- 
face. 

From these considerations it appears that the conclusion of refer- 
ence 10 that inviscid flow along streamlines downstream of the nose of 

a It is interesting to note that in ideal gas flows, r becomes just a 
function of these variables as the value of the hypersonic similarity 
parameter becomes large compared to 1 (see work of Oswatitsch, refer- 
ence 33 ) noting that his results can readily be extended to three- 
dimensional ideal gas flows using the characteristics method of this 
paper). 
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noninclined bodies of revolution traveling at high supersonic speeds is 
often of the Prandtl-Meyer type (in regions free of shock waves) may 
apply also to other steady three-dimensional flows. It is true, too, 
that in the latter case, just as in the former case, this conclusion is 
consistent with the predictions of the hypersonic similarity law for 
steady flow about slender shapes . 

One question remains to be considered, namely, where do the stream- 
lines go in the disturbed flow? To clarify this matter, it is conven- 
ient to study further the implications of equation (25). For this purpose 
we combine equation (25) with the transformation equation 


8A _ M f 8a 8a ''N 

8y 2 \5Ciy dCay/ 


to obtain the relation 
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From this relation we deduce either that to the order of a number (curva- 


ture) small compared to 
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Equation (27) implies vortical flow, however, which type of flow cannot 
be treated by the present analysis since equation (25) is violated. 9 
Equation (28) is then the requirement consistent with the basic assumptions 


9 This conclusion Is particularly evident in the case of pure vortical flow, 
or say vortical flow with a superimposed uniform stream directed along 
the axis of the vortex, in which cases 8p/8x = 0, and hence equation (27) 
certainly does not follow from equation (22). 
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of this analysis. Comparing the relations of equation (28) with the 
transformation equation 

dA _ M / 8A \ 

^ 2.J¥^L V Sc iy + Sc 2y/ 

leads one to the conclusion, however, that 


8s 

» 


8x 


8x 


or, in effect, that surface streamlines can, to the accuracy of this 
analysis, be considered geodesic lines of the surface. With this Inf or - 
tion we are enabled to construct the flow field about a body, having 
once determined,, for example, the flow in the region of the leading edge 
(or edges) thereof. This result follows since a geodesic line, and hence 
a streamline, on the surface is fixed, provided its direction at any 
point is given (see, e.g., reference 3^). 10 With this knowledge of the 
location of surface streamlines, flow in the planes tangent thereto and 
normal to the surface may be calculated approximately in the relatively 
thin region between the surface and bounding shock waves, using the gen- 
eralized shock-expansion method after the manner described in reference J. 

* A partial check on these observations is afforded by studying the 
flow about a swept airfoil. In this case flow at the surface may be 
calculated with good accuracy, using the shock-expansion method in com- 
bination with simple-sweep theory. For thin airfoils (on the surface of 
which the appropriate geodesics have essentially the direction of the 
free stream) the extended shock-expansion method of this paper reduces to 
the slender-airfoil method of reference 7* Thus, in this case, it is 
evident from the results of reference J that the extended method will 
predict surface pressure coefficients in error by less than 10 percent, 
providing the component of free-stream Mach number normal to the leading 
edge is greater than about 3« It, is of interest also to consider a thick 
airfoil to ascertain the accuracy with which this method applies to flow 
with appreciable curvature. To this end, surface pressure coefficients 
and streamlines have been calculated for a 20-percent-thick biconvex 
airfoil (at zero incidence) swept 6o° and operating at Mach numbers of 10 
and infinity (7 = 1.4). The results of these calculations are presented 
in figure 1, and it is observed that the pressure distributions determined 
with the shock-expansion method for swept airfoils and the extended shock- 
expansion method are in reasonably good agreement at both Mach numbers. 


10 If a sudden change of surface slope causes an oblique shock wave or a 
concentrated Prandtl-Meyer type expansion fan, the stre aml ines in the 
downstream direction are defined on the basis of their flow direction 
immediately following the discontinuity in slope. 
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The streamlines are also in reasonably good agreement over the forward 
portion of the airfoil, although, as would be expected, somewhat poorer 
results are obtained over the afterportion. It is not surprising, in 
view of the underlying assumptions of the extended shock- expans ion method, 
that it is generally more accurate at the highest Mach number. 


In the preceding discussion circumstances were deduced under which 
steady flow at high supersonic speeds about three-dimensional shapes 
could be constructed approximately, using the basic tools of two- 
dimensional supersonic flow analysis, namely, the oblique shock equations 
and Frandtl-Meyer equations. Several possible exceptions to these cir- 
cumstances immediately come to mind. These include coni cal -type flow 
and flow in the region of the tip of a wing, or at the discontinuous 
juncture of a wing and body, to mention a few. In such flows equation (25) 
may not be satisfied, in which case two-dimensional flow in planes normal 
to a surface cannot be expected. 11 It might be reasoned, therefore, that 
these flows cannot, in general, be treated by the proposed method. This 
observation may be correct; however, in the one case investigated thus 
far in this connection, namely, flow in the region of the nose of non- 
lifting bodies of revolution (see reference 10), it was found that 
although equation (25) is not satisfied, flow along streamlines is never- 
theless of approximately the Prandtl-Meyer type. Thus we are led to 
expect that perhaps a less restrictive requirement than the satisfying 
of equation (25) may be imposed to insure that flow along strea m lines is 
of this type. Such a requirement is in fact easily obtained by reconsid- 
ering equation (22) in the form 


dp _ pU 2 dS M f dS ^ 1 dA ^ 

d* ~ Vm 2 -i L ” </m 2 -i V Sc iz m dy/ _ 

thus yielding 


dS 


db 

1 dA 

dx 

VM 2 -! 

^ c lz 

| - 
M dy 


(30) 


(31) 


It is evident that equation (3l) embraces equation (25) as a special case 
and that Prandtl-Meyer flow obtains along streamlines if 


d5 ~ 1 dA 

dC lz M dy 


(32) 


■^One may note that in seme cases of this nature, the flow in osculating 
planes of the streamlines may be of the two-dimensional or even the 

simpler Prandtl-Meyer type, although these planes may not be normal to 
the surface . 
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to the order of a number small compared to 


«/m 2 -i 

M 


85 

8x 


This result 


implies that although flow inclination angles are not necessarily con- 


stant along 
tion (17)). 


'lz 


lines, pressure is approximately constant (see equa- 


It is clear that the increased generality of the above result has 
been obtained at some expense in our knowledge of the streamline flow 
pattern. For example, it is not now indicated that (within the framework 
of this analysis) surface streamlines may generally be taken as geo- 
desics - additional knowledge of the flow must be had in order to deter- 
mine these streamlines. If they are known, however, the calculation of 
the whole flow field is materially facilitated by the above considerations 
To illu strate, consider a nonlifting body of revolution (see sketch) 



for which we assume that flow at the vertex is known either from say 
reference 8 or reference 10. The meridian curve of the body is broken 
up into short segments as shown, and flow is constructed along the first- 
family Mach lines emanating from the intersection of these segments. The 
requirement to be satisfied is that the pressure change across these lines 
be c ons tant along their length. 12 The construction proceeds then in a 
manner analogous to that for the two- dime nsional airfoil discussed in 
reference 7* 

Thus far only steady flows have been considered. The problem 
naturally arises of extending these considerations to nonsteady flows. 

Some aspects of this matter will now be discussed. 


Nonsteady Flow 


The methods of analysis in this case are entirely analogous to those 
employed in the study of steady flow, the singular contrasting feature 


1 p 

In this manner small changes in pressure along lines can be 

accounted for approximately in the predominantly conical flow near the 
vertex of the body. 
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being that derivatives with respect to time in equations (l) through (ll) 
cannot nov be neglected. With this point in mind, only pertinent results 
are discussed below. 


Characteristics method .- The compatibility equations relating fluid 
properties along Mach lines may be written as follows: 



85 ] 1 /c&\ + 1 
8c lz M Vdy / U 


V 


« 

Vm 2 -1 /85\ + JL /8u\1 

M at ) pu 2 V dt J ” MU \ 8 t J J 


( 33 ) 

and 


8p _ + P U g 

Sc 2z a/m 2 -1 


85 1 f 8a\ + l 

8c 2Z m \8yy U 


//m 2 -! ( 8s\ M ( 8p\ • 1 / 8u V 
m \ 8ty pu 2 \ 8ty mu v Sty_ 

( 34 ) 


The definition of the x - z plane as the osculating plane of a pathline 
(streamline in steady flow) remains as before, hence equation ( 19 ) still 
applies in the x - y plane in the region of the origin. The rotation 
of the osculating plane and variation of the principal curvature of a 
pathline with motion along it are now, however, obtained with the aid of 
the relations 


( 35 ) 

dx 2 pU 2 dx \5y / 3z \ dx / K 

and 

22 „ - _i_ A f h') - f2 + 32E + JLs!E N ) 2 ( 36 ) 

dx 2 pU 2 dx \8z/ \ 8z U dx pa 2 dx/ dx 

where 

sri + !Gl) (37) 

These equations in combination with the energy and state equations are 
employed in the same manner as in the. case of steady flow to construct a 

flow field, proceeding from an initial value surface. It is clear, how- 

ever, that because of the unsteady nature of the flow, this surface is 
not necessarily fixed in space, nor are fluid properties necessarily 
constant on it. Thus, in order to construct the flow, it will, in general. 
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be necessary to start new osculating planes from this surface at short 
intervals of time, each plane being attached to a particular element of 
fluid as it moves through the field. By way of comparison, then, we 
recall that in going from steady two-dimensional to three-dimensional 
flow with the characteristics method, it was necessary to construct the 
flow in a family of surfaces (located adjacently in say the spanwise 
direction) rather than just a single surface. Analogously, in going from 
three-dimensional steady to three-dimensional nonBteady flow, it is neces- 
sary to construct the flow in a family of spaces located adjacently in the 
"time" direction, rather than in just one space. Quite obviously such a 
series of calculations poses so formidable and time consuming a problem 
as to be questionably feasible at present; hence they will be considered 
in no greater detail here. Rather, let us turn our attention to the 
approximate method of calculating nonsteady flows. 

Approximate method .- As in the case of the corresponding steady-flow 
analysis, it is convenient to consider the expression for pressure gra- 
dient along a pathline. Thus we have 

( i-Dz ^ 1 J 5 a 1 [ 1 f Sp^ 8u 

dx uW J \1+D Z / 1 8y U L pU 2 Veit/ St 

( 38 ) 

where now the x - z planes are taken normal to the surface swept out by 
elements of fluid moving along the body. Upon inspection of this relation 
and the compatibility equations, it becomes clear that the critical 
requirement for two-dimensional flow of the generalized Prandtl-Meyer type 
in these planes is, in addition to the one previously derived from steady 
flow considerations, that 

1 85 <<; dS 

U St dx 

or 

_1_ S6_ d5 

M° Sxq dx 

where 


M° = — . Xq = aot 
9-0 

a Q being the velocity of sound in the undisturbed stream. Now so long 
as the Mach number of the undisturbed stream and the local Mach numbers 
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are large compared to 1 , M° is large compared to 1 since the speed of ... 
the undisturbed stream and the speed of the local flow cannot differ 
greatly. Thus, with this restriction, the requirement expressed by 
equation ( 39 ) is simply that the induced curvature of the flow | 8 &/ 8 x 0 | , 
associated with the nonsteady motions of the fluid, cannot exceed in 
order of magnitude the total curvature |d5/dx| of the flow. Provided 
this is the case, and equation ( 25 ) is satisfied, equation ( 38 ) reduces, 
of course, to 



One observes that this result is not applicable to as wide a class of 
shapes as that for steady flow, since the local Mach number of the dis- 
turbed flow is now required to be everywhere large compared to l . 13 
This additional requirement manifests itself, since otherwise, nonsteady 
disturbances created an appreciable distance upstream and/or downstream 
of a particle could significantly influence its behavior in the disturbed 
flow field (see sketch, noting that in case of thick body, particle b 
is influenced by disturbances originating in particleB a and c). 



Thus equation (4o) applies only in disturbed flow fields about thin or 
slender shapes (i.e., shapes producing flow deflection angles small com- 
pared to l). In such cases, pathlines in the surfaces swept out by 
elements of fluid adjacent to the shapes are approximated by geodesics 
or, even simpler, lines of curvature of these surfaces. It is not to be 
implied, of course, that pathlines must always be such curves in order 
for fluid properties to behave as in Prandtl-Meyer flow. In fact, just 
as In the case of steady flow. If equation (3l) rather than equation (25) 
is satisfied, pathlines are not necessarily geodesics (or lines of curva- 
ture) although equation (39) and hence equation (40) hold along these lines. 

13 The net simplification of requiring only that the hypersonic similarity 
parameter of the flow be large compared .to 1, Is, in general, that flow 
in osculating planes may be treated as nonsteady, two-dimensional. 
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One notes that within the framework of this approximate analysis, 
the calculation of nonsteady flows at least at the surface of slender 
bodies traveling at high supersonic speeds should not prove unduly diffi- 
cult. To illustrate, consider an oscillating airfoil as shown in the 
sketch: 


Shock wav e 


Direction of rotation 
of airfoil 


The pressure at any point along the pathline shown is readily deduced by 
simply integrating equation (40) along this line from the leading edge 
of the airfoil to the point in question. The whole flow field as a func- 
tion of time may be calculated by employing the generalized shock- 
expansion method for steady flows (see reference 7) in a series of planes 
located small distances apart in time. This example serves to emphasize 
that the time history of fluid elements must be known, at least to the 
extent of fixing their initial flow direction and entropy. It is also 
evident that again, as in the case of steady flow, the general results 
of the analysis are consistent with the predictions of the hypersonic 
similarity law for nonsteady flows about slender related shapes (refer- 
ence 5) . 


Airfoil at time t2 




Pathline of particle striking 
leading edge at time fj 

Airfoil at time t/ 



CONCLUDING REMARKS 


A method of characteristics for solving steady three-dimensional 
supersonic flow problems has been considered. It was found that compat- 
ibility equations relating fluid properties along characteristic lines 
could be obtained in a simple form by employing pressure and flow incli- 
nation angles as dependent variables. No significant restrictions were 
imposed on either the equation of state obeyed by the fluid, or the rela- 
tions defining its specific heats. These features of generality were 
retained for the specific purpose of enabling more accurate application 
of the method to the calculation of flow fields about missiles traveling 
at high supersonic airspeeds . Such application requires , of course, a 
predetermined knowledge of fluid properties along some surface in the 
disturbed flow. Extension of the method to treat nonsteady flows was 
considered briefly. 

It was also undertaken to obtain an approximate method for calculat- 
ing flows about bodies traveling at high supersonic speeds. It was found 
that when the flight Mach number is sufficiently large compared to 1, flow 
in the osculating planes of streamlines in regions free of shock waves 
may frequently be of the generalized Prandtl-Meyer type - surface 
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streamlines in this event may be taken .as geodesics. In the case of 
slender shapes, these results apply to nonsteady as well as steady flows, 
provided the induced curvature of streamlines does not exceed the total 
curvature in order of magnitude. It is concluded from these and other 
considerations that two- dime hsional-flow equations may be applicable to 
a relatively wide class of flows, and hence configurations, at high 
supersonic speeds . 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Aug. 15, 1952 
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APPENDIX 

SYMBOLS 


a 

c 

Qlzj C 2Z 


Cp 

M 


local speed of sound 

chord of airfoil (measured normal to leading edge) 

characteristic coordinates in x - z plane (Ci z is positively 
inclined with respect to x) 

pressure coefficient ( — B ~ E 0 — ) 

VpoUo 2 /2/ 

Mach number (ratio of local velocity to local speed of sound) 


M° Mach number (ratio of local velocity to speed of sound in the 

undisturbed stream) 

p static pressure 

s entropy 

t time 


U,V,W components of fluid velocity along the x, y, and z axes, res- 
pectively 

x,y,z rectangular - coordinates 

7 • ratio of specific heat at constant pressure to specific heat 

at constant volume 


6 angle between x axis and tangent to projection of streamline 

(or pathline) in x - z plane 

A angle between x axis and tangent to projection of streamline 

(or pathline) in x - y plane 

p density 


Subscript 


o free-stream conditions 
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Figure /.- Comparison of surface streamlines and pressure distributations calculated with the extended shock -expansion 
method and the 
angle = 60 °). 
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